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Abstract 

One of the simplest examples of a T^T-symmetric quantum system is the 
scaling Yang-Lee model, a quantum field theory with cubic interaction and 
purely imaginary coupling. We give a historical review of some facts about 
this model in d < 2 dimensions, from its original definition in connection with 
phase transitions in the Ising model and its relevance to polymer physics, to 
the role it has played in studies of integrable quantum field theory and of VT- 
symmetric quantum mechanics. We also discuss some more general results 
on T^T-symmetric quantum mechanics and the ODE/IM correspondence, 
and mention applications to magnetic systems and cold atom physics. 



1 Introduction 



The scaling Yang-Lee model in two dimensions is one of the simplest examples of 
an interacting integrable quantum field theory. It has often been used by the inte- 
grable models community as a testing ground for new mathematical and numerical 
tools. Aspects of the theory of exact S-matrices [HIS]) the thermodynamic Bethe 
Ansatz (TBA) [3], the truncated conformal space approximation (TCSA) [11[5], 
the form factor approach for correlation functions ^ X\ , the exact off-critical g- 
function [U [9] , the excited state and the boundary TBA methods [TOl ttH El [5] 
are all techniques which were very successfully applied in their early stages to the 
study of non-perturbative phenomena in this simple model. 

Deep in the ultraviolet regime the properties of the scaling Yang-Lee model 
are governed by the conformal field theory A^2,5 [121 [T3] which has a negative 
central charge and a single relevant spin-zero field (p, with negative conformal 
dimensions. The off-critical integrable version of the model corresponds to the 
perturbation of the conformally-invariant A^2,5 action by the operator (p with a 
purely imaginary coupling constant [l]. 

Despite its apparent lack of unitarity, the Yang-Lee model in d space-time 
dimensions is at least relevant in condensed matter physics: it is related to the 
theory of non-intersecting branched polymers in d -|- 2 dimensions ^4J . 

Considerations of the Yang-Lee model led Bessis and Zinn- Justin to conjec- 
ture the reality of the spectrum of the Schrodinger equation with cubic potential 
and purely-imaginary coupling, which in turn prompted Bender and his collab- 
orators to a more general study of Schrodinger problems with complex poten- 
tials \15\ I16j . The paper [15] marks the beginning of "PT-symmetric quantum 
mechanics and quantum field theory as an area of intensive research. This short 
review and the corresponding introductory lecture given in Mumbai are especially 
addressed to students approaching the field of quantum mechanics in the com- 
plex domain. However, the brief discussions on conformal field theory, condensed 
matter physics, integrable models and their correspondence with the theory of 
ordinary differential equations may also stimulate more experienced researchers 
to look into some of the many open problems related to non-unitary integrable 
quantum field theories. 

2 The Yang-Lee edge singularity 

Consider the ferromagnetic Ising model with non-zero external magnetic field H 
n[s,i{] = -Y,SiSj+ii^Si, Si = ±l (2.1) 

<ij> i 
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where the sum is over nearest-neighbour sites in a 2d square lattice. The partition 
function is 



In the thermodynamic limit, at H = with T above a certain critical temperature 
Tc, the system is in a disordered phase with zero magnetisation M(r, H)|h=o = 
(si) = 0. At T < Tc the system is instead in an ordered phase where the Z2 
symmetry is spontaneously broken and M(T, H)|h=o 0- Finally at H = and 
T = Tc the model undergoes a second order phase transition where typical spin 
configurations exhibit fluctuations at all length scales and the continuum limit 
version of the model is invariant under conformal transformations. 

A result that dates back to Yang and Lee in 1952 [n\ 118] states that the 
partition function Z{T, H) of the finite lattice version of the Ising model, at fixed 
T, is an entire function of H with all the zeros located on the purely imaginary 
H-axis. In the thermodynamic limit the zeros become dense, covering the entire 
imaginary axis from —00 to 00 apart from a possible single finite gap centered 
at H = 0. When T < Tc, there is no gap in the zero distribution and the 
magnetization M(T, H) has a finite discontinuity as H crosses the origin along the 
real axis. Correspondingly the system undergoes a first-order phase transition. 

For T > Tc there is a gap for |II| < hc{T). The edges H = zizihc{T) of this 
gap are branch points for the magnetization M(r, H) and according to Fisher [19] 
they can be considered as conventional critical points. In higher dimensions they 
correspond to the infrared behavior of a scalar field theory with action 



The action (j2.3p is related to the Parisi-Sourlas theory of non-intersecting branched 
polymers in D = d + 2 dimensions [14] . 

3 The Yang-Lee model and branched polymers 

Following [141 1201 [2T] . non-intersecting branched polymers in D dimensions cor- 
respond to the classical field equation 




(2.2) 



Si=±l 




(2.3) 



(3.1) 



where the field potential V{(j)) is 



v{<f>) = -{h-hc)m + ^4>''ix) 



(3.2) 



and ^(x) is a stochastic (white noise) variable with 



{amy)) 



{x - y). 



(3.3) 
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Figure 1: A branched-polymer tree diagram. 

The polymers can be visualized using tree Feynman diagrams (see figured]) with 
corresponding Feynman rules: 

• A segment in the graph corresponds to the Green's function of — . 

• Each vertex carries a factor g . 

• Each free end carries a factor —{h — he) + 
The n-point correlation functions are defined as 

{cPixi) . . . ^{xn)) = Jv^ 4>{xi) . ■ ■ ^{xn)e--2f^'°^"'^''> , (3.4) 

where the 4'{xi)^s solve the classical field equation (|3.1|) . 

Resolving the constraints using the Faddeev-Popov method introduces a Grass- 
mann ghost field and reveals that the theory has an unexpected invariance under 
supersymmetric transformations. A consequence of this hidden supersymmetry 
is a surprising dimensional reduction property: correlation functions whose ar- 
guments are restricted to a d = D — 2 subspace are the same as those for the 
Yang-Lee theory in d dimensions. 

The simplest case corresponds to the branched polymer problem in D = 2 
dimensions. In d = D — 2 = the gradient term in the Yang-Lee action (12. 3p 
is absent and the vacuum-vacuum functional integral becomes, after a complex 
rotation, 

|^A'(^)=£ 

This is one of the very few exactly-known scaling functions in = 2. This 
model describes the crossover from self-avoiding loops (vesicles) to branched poly- 
mers [22\ [20| |2T| where p ~ is the pressure difference between the inside and 
the outside and (xc — x) ~ (/i — hc)g~^ is related to the monomer fugacity in 
the discretized version of the model. This result can be interpreted from another 
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Figure 2: The self-avoiding loops to branched polymers crossover transition. 



perspective: since the original polymer problem is defined in D = 2 dimensions 
and the Airy function is the solution of a Id Schrodinger equation with linear 
potential, we have a beautiful physically- motivated confirmation of an ODE/IM 
correspondence p3 i linking the spectral theory of simple ordinary differential 
equations and the theory of integrable models. 

4 Conformal field theory in two dimensions: the Ising 
and the Yang-Lee models 

From section [2] we know that the Ising model at H = and T = Tc is critical and 
invariant under conformal transformations. The edge points H = zizihc{T) can 
also be considered as conventional critical points ^19j where the model exhibits 
conformal invariance [12]. From a renormalisation group perspective the point 
(H = 0, T = Tc) and the pair of equivalent points (H = ±ihc{T), T > Tc) 
correspond to two distinct universality classes. 

Two dimensional conformal transformations coincide with analytic transfor- 
mations 

z^fiz) , z^m (4.1) 
with z = X + iy, z = x — iy. The corresponding infinitesimal generators are 

l^ = -z^+'d, , [„ = -z"+ia,- , (4.2) 

and satisfy the Witt algebra 

[ImJn] = (m- n)lm+n , [Imjn] = {jTl - n)lm+n , [L, ^m] = . (4.3) 

At the quantum level the geometry is important and commutation relations in- 
volving infinitesimal generators that change it become anomalous. The Witt alge- 
bra is replaced by the Virasoro algebra Im Lm , Im — *■ Lm, where [L„,Lm] = 
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and 

[Lm, Ln] = (m - n)Lm+n + ]^("'^ ~ n)5n+mfl , 

[Lm-, Ln] = {m- n)Lm+n + ]^("'^ ~ n)5n+mfl , (4.4) 

with c the Casimir coefficient describing the couphng of the quantum fluctuations 
to the finite geometry of the system. The quantity c is also cahed the central 
charge, or conformal anomaly. 

4.1 The Kac table 

In two dimensions a universality class is completely identified by its central charge 
and a complete set of operators, or fields, with their operator product expansions. 
Among all the operators the primary operators play a central role: they are the 
ones that transform in the following simple way under a conformal transformation 

<t^{zrz) - {d,f{z))\dJ{z)fct>{f{z)J{z)) , (4.5) 

where h and h are the holomorphic and the anti-holomorphic conformal dimen- 
sions respectivelyB. The remaining operators - called descendants - can be ob- 
tained from the primary operators by acting on them with combinations of L—n 
and L-n. For generic values of c there are infinitely-many primary operators, but 
in their pioneering paper [21] on conformal field theory, Belavin, Polyakov and 
Zamolodchikov showed that for each pair of coprime integers p < q, a 'minimal' 
conformal field theory A1p,q can be defined, with 

c = l-6(P^<l (4.6) 

pq 

and a finite number of primary operators, with conformal dimensions 

K,.= ^-^^^^^:l^i^, (l<.<p, l<.<q./p). (4.7) 
4pq 

The set of conformal dimensions forms the so-called Kac table of the model. The 
critical Ising model corresponds to the minimal model ^^3,4 with central charge 
c = 1/2 and Kac table 

1 <-> hi^i = hi^i = , e <-> /ii,3 = /ii,3 = 1/2 , a hi^2 = hi^2 = 1/16 , (4.8) 

where 1 is the identity operator, e and a are respectively the energy and the mag- 
netic field operators, while the Yang-Lee edge singularity corresponds to M.2,5 , 
with central charge c = —22/5 and Kac table 

I ^ hi^i = hi^i = , (/) ^ /ii^2 = /ii,2 = -1/5 . (4.9) 

Notice that the second of these conformal dimensions is negative, a refiection of 
the non-unitary nature of the Yang-Lee model. 



* In the following we are only interested in spin-zero operators with h — h. 
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5 The Scaling Yang-Lee model in two dimensions 



A 2d conformal field theory is a particular example of an integrable quantum field 
theory, the integrability of the model being a direct consequence of the existence 
of the infinite number of commuting conserved charges built from {L„,L„}. It 
was an idea of Sasha Zamolodchikov [25] that a conformal field theory can be 
perturbed or driven away from the critical point by one of the relevant spin-zero 
fields of the theory, in such a way that integrability would be preserved. For 
particular primary fields he was able to deduce that an infinite set of commuting 
conserved charges survives after the perturbation and that in these cases the 
model remains integrable. It is now known that both (pi2 and (pi^ perturbations 
of the minimal models lead to integrable off-critical quantum field theories. The 
perturbation of A^2,5 by the (p operator 



is therefore integrable, and can be studied using powerful techniques from the 
theory of integrable models. Results that have been obtained include: 

• The exact 2-body S-matrix was proposed by Cardy and Mussardo [Ij and 
independently by Smirnov [2j in 1989. The resulting theory for (h — hc) > 
contains a single particle species with a cp^ interaction. No other bound 
states exist. 

• Based on the S-matrix of [Tl[2], Aliosha Zamolodchikov derived a set of ther- 
modynamic Bethe Ansatz (TBA) equations describing the exact ground- 
state energy of the model on an cylinder with finite circumference R [3]. 

• In 1996 the thermodynamic Bethe Ansatz equations were generalised to the 
full set of excited states [101 E] • 

• In 1990 Aliosha Zamolodchikov found an exact form-factor expansion for 
particular two-point correlation functions [6]. 

Finally, the low-lying energy levels of the theory can be studied numerically 
using another method, proposed by Yurov and Aliosha Zamolodchikov [1], called 
the truncated conformal space approximation. This relies on the diagonalisation 
of the perturbed Hamiltonian defined on a circle of circumference i? on a suitable 
truncation of its Hilbert space, using the basis provided by the (exactly-known) 
space of states of the unperturbed conformal field theory. Figures [3] and U] repro- 
duce some of the results found by Yurov and Zamolodchikov in 1990. Figure [3] 
shows the spectrum of the Yang-Lee model at fixed (h — hc) > 0, as the circumfer- 
ence R of the cylinder is increased. One can see that the spectrum is completely 
real. The bottom line represents the ground-state energy Eq and above this line 
lies the single one-particle state Ei, the energy gap tending asymptotically to the 




(5.1) 
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mass m of a single neutral particle at rest. Above this second line, asymptotically 
again at E2 — Ei = m there is the lowest two-particle excited state, and above it 
the continuum. 



Re E(R) 




Figure 3: Yang-Lee energy levels for {h — he) > 0. 

The situation for (/i — he) < is very different. For small values of R all 
the energy levels are still real, but the ground and the first excited states meet 
at a critical value Rc , beyond which they form a complex-conjugate pair. At 
the special point where this happens, there is a square-root branch point of the 
ground-state energy. 
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Re F{R) 




Figure 4: Yang-Lee scaling functions F{B) = RE{R)/tt for {h— he) > (left-hand plot) 
and {h — he) > (right-hand plot). 

Figure m shows this and other mergings, and illustrates the contrast between 
the situations for the two signs of {h — he)- The plots show the so-called scaling 
functions F(R) = RE{R)/tt. Conformal field theory gives the (finite) limits of 
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these functions as R ^ 0. For R ^ the scahng functions move away from these 
limits (in fact, as analytic functions of {h — hc)R^'^^^), and for (h — hc) < certain 
pairs meet at square-root branch points, in a way that will see a striking echo in 
a quantum-mechanical context in the next section. 

As a more technical remark, note that the explicit factor of i appearing in the 
action (jS.ip is a result of the particular normalisation of the perturbing field (j){x) 
used in [4j , which ensured that the corresponding primary state in the conformal 
field theory had positive norm-squared. However there is no compelling reason 
to make such a choice: if Ln = L-n, then it follows from the Virasoro algebra 
(|4.4p that, if 1-0) is the state corresponding to a primary field tp{x), then the 
squared norm of, for example, L-n\'4') is 2n(/i+ -^{ri^ — ^)) times that of {(p). In a 
non-unitary conformal field theory such as the Yang-Lee model, where both c and 
some of the conformal dimensions h are negative, there is therefore no natural 
way for all (primary and descendant) states to have positive squared norms, and 
from many points of view it is preferable to modify the normalisations chosen 
in [4j, so that all coefficients in operator product expansions are real (see, for 
example, [26l [5l [27] )• This removes the factor of i from ()5.ip . but of course the 
theory remains non-Hermitian, since the norm in the underlying Hilbert space 
is not positive-definite. A similarly non-positive-definite norm will appear in a 
quantum-mechanical context in the next section. 

Finally for this section, we should mention that VT symmetry can also be 
found in integrable lattice models, where the continuous dimensions of space and 
time are replaced by a discrete lattice. An initial set of examples was discussed 
in [28], but it is clear that many further cases remain to be studied. 

6 PT-symmetric Quantum Mechanics 
6.1 Spectral reality 

In the early 1990s, motivated by the Yang-Lee edge singularity, Bessis and Zinn- 
Justin considered the spectrum of the non-Hermitian Hamiltonian 

HBZj=P^ + ix^ (6.1) 
defined on the full real line. The corresponding Schrodinger equation is 

- -r^il'ix) + ix^ij{x) = Eni^ix) , ip{x) € L'^{1R). (6.2) 

Perturbative and numerical work led them to conjecture [29] 
• The spectrum of Hbzj is real, and positivell] 

^ While the Bessis - Zinn- Justin conjecture initiated the line of work we want to describe 
here, similar observations had in fact been made before. For example, in 1980 the oscillators 
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In 1997, Bender and Boettcher [T5] proposed an apparently-simple generalization: 



Hm =P^ - {ixf^ , M G R , M > . (6.3) 

To render the 'potential' — (ix)^*^ single- valued for all values of M, a branch cut 
should be placed along the positive imaginary axis of the complex x plane. For 
M < 2, the analytic continuation of the spectrum at M = 3/2 (Bessis and Zinn- 
Justin's original problem) can be recovered by imposing the boundary condition 
'ip{x) £ L^(]R), or equivalently demanding that ■ip{x) should tend exponentially 
to zero as \x\ — > oo along the positive and negative real axes. However, when M 
reaches 2, the potential becomes — and the eigenvalue problem as just stated 
changes its nature dramatically. To see why, consider the WKB approximation to 
the wavefunction as x tends to infinity along a general ray in the complex plane: 



V'i ~ P ^^^"^ 6Xp 



.^je{i+M) 1+M 

■ A/+1 ^ P 



(x = pe'^/i) . (6.4) 



For most values of 6, one of these solutions will be exponentially growing as p ^ oo 
(a dominant solution) and the other exponentially decaying (or subdominant). 
But when 6 is such that 



neither solution is dominant: both decay algebraically. The angles 
„ vr 3-7r Svr 



"2M+2 ' 2M+2 ' 2M+2 



(6.6) 



define anti-Stokes lines, which divide the complex plane into Stokes sectors. 
Across an anti-Stokes line, the dominant and subdominant solutions swap roles. 

At M = 2, the positive and negative real axes coincide with anti-Stokes lines, 
and the original eigenproblem ceases to have a discrete spectrum. The issue 
can be resolved if the perspective is widened by allowing the wavefunction to be 
defined along general contours in the complex plane, rather than just the real 
axis. Two types of eigenvalue problem are natural. Lateral problems correspond 
to boundary conditions which require wavefunctions to decay at large |x| along 
a contour that begins and ends in a pair of Stokes sectors. By contrast, radial 



p + X + igx were studied ^0], and the reality of individual eigenvalues proved for g real 
and \g\ sufficiently small. However it is important to note that this does not suffice to prove 
that the whole spectrum is real for any non-zero value of g. The reason is that the 'sufficiently 
small' value of \g\ may depend on the particular eigenvalue under consideration. Without a 
global bound away from zero for these values, reality can still be lost sufficiently high in the 
spectrum for all non-zero values of g, even if each individual eigenvalue ultimately becomes real 
as \g\ — > 0. (This is exactly the behaviour seen in Bender and Boettcher's generalisation of the 
Bessis - Zinn- Justin problem as M — > 1~, illustrated in figure [7] below, and it is not ruled out 
by the arguments of [3D] ) . A full proof of spectral reality for -\-x'^ -\- igx^ - in fact valid for any 
real value of the parameter g - was finally provided by Shin in 2002 [31] . following the approach 
of 132] . 
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X 

Im ^ 




Figure 5: Sample quantisation contours. The dashed hues separate Stokes sectors. 

problems correspond to a quantisation contour which starts at x = 0, and runs out 
to infinity in a given Stokes sector. Figure [5] depicts some sample quantisation 
contours. Some eigenproblems defined in this way will be related via simple 
changes of variables, but those that are not have completely different spectra. 

For M < 2, the Bender-Boettcher problem (16. 3p is an instance of the lateral 
problem with the (next-nearest-neighbour) Stokes sectors that include the real 
axis. For M > 2, these sectors rotate down in the complex plane as shown in 
figure O and no longer include the real axis. Demanding that the quantisation 
contour begins and ends inside this pair of Stokes sectors for all positive values 
of M gives the analytic continuation of the Bender-Boettcher problems as M 
increases past 2. 

Im I X 

A 




Figure 6: The Bender-Boettcher contour, for M just larger than 2. 
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Numerical and analytical evidence led Bender and Boettcher to conjecture 
that, with these boundary conditions understood, the spectrum of (j6.3p is real 
and positive for all M > 1, despite the non-hermiticity of the problem. The real 
eigenvalues as a function of M are shown in figure[71 As soon as M decreases below 
1, infinitely-many eigenvalues pair off and become complex, while the spectrum 
at M = 1/2 is empty. 




0.5 1 1.5 2 2.5 3 
M 



Figure 7: Real eigenvalues of — (ix)^*^ as a function of M. 

For M = 2, the upside-down quartic potential along the complex contour 
described above, the reality of the spectrum can be understood via a spectrally- 
equivalent Hermitian Hamiltonian [33^ [34l [35] . Jones and Mateo found a partic- 
ularly transparent way to derive this result [35] using the simple variable change 
X = —2iy/l + iw, which maps the complex contour onto the full real line. This 
variable change followed by a Fourier transform turns the "PT-symmetric Hamil- 
tonian H2 defined on the contour in figure [6] into 

H2=p^ + - 2w (6.7) 

defined on the full real line, which is explicitly Hermitian. Apart from M = 2 
and M = 3 02], mappings to Hermitian Hamiltonians have not been found for 
the Bender-Boettcher problems Hm for other values of M > 1. 

Bender and Boettcher interpreted the phase transition from entirely real eigen- 
values to infinitely-many complex eigenvalues as a spontaneous breaking of VT 
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Figure 8: Real eigenvalues of -ffi — + x^^ (ixY as a function of e (plots taken 
from [IS])- 

symmetry. The parity reflection V and time-reversal T operators act as follows: 

V : X ^ —X , p —p 

T : X ^ X , p ^ —p , i — > —i. 

The combined operator VT commutes with the Hamiltonian Hm for all values 
of M. The eigenvalues of VT are pure phases and with a multiplicative rescaling 
of the eigenfunctions they can be set equal to 1. If the eigenfunction tpn of Hm 
is also an eigenfunction of VT, then we say that the PT-symmetry is unbroken. 
In this case, using commutativity we have 

VT Hm^u = VT Enipn = K 'PT^n = E* (6.8) 

Hm VTi^n = Hm = ^„^n (6.9) 

and so En is real. Unfortunately, this does not seem to offer a quick route to a 
proof of spectral reality, since it is a non-trivial task to prove that VT symmetry 
is unbroken for all wavefunctions. The only proof currently known of the reality 
of the spectrum of Hm for general M > 1 [32] was obtained by a rather different 
route, using ideas from the ODE/IM correspondence [23]. 

Spectra for problems involving next-next-nearest-neighbour Stokes sectors and 
next-next-next-nearest-neighbour Stokes sectors are shown in figure [HI These 
plots, taken from [16], show the real eigenvalues as a function of e for the K = 2 
and K = 3 cases of 

H^/) = p2 + x'^^iixf {Ken, e real) , (6.10) 

with the implicit understanding that in each case the quantisation contour at 
e = is the real line. The spectrum aX K = 2, shown in figure [S^i, is conjectured 
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to be real for all e >0, and undergoes an extra 'phase transition' at e = —1. The 
real eigenvalues for next-next-next-neighbour Stokes sectors (the K = 3 case of 
()6.10p ) are shown in figure [Hb- This time there are phase transitions at e = — 2, — 1 
and 0, and the spectrum is conjectured to be real for all e > 0. For integer values 
of e > 1 — iiT, the reality of the spectrum of (|6.10p was proved in [36], using 
techniques related to those of [32J. (The negative values of e covered by the proof 
correspond to the extra phase transitions just mentioned.) 

An alternative to looking at eigenproblems in varying pairs of Stokes sectors is 
to stay with the same Stokes sectors as for the original Bender-Boettcher problem 
and explore other possibilities for the potential. In [37] the effect of the addition 
of an angular-momentum-like term l{l+l)/x'^ was investigated, the Hamiltonian 
becoming 

Hm,i =p'- {ixf^' + l{l + l)/x^ ■ (6.11) 

The results are shown in figure [3 note the completely reversed connectivity of 
the real levels in figure [9^, compared to that of figure [71 The remaining figures 
show how this connectivity changes as / increases toward zero, driven by the lowest 
eigenvalue. This behaviour was subsequently understood analytically, again using 
tools borrowed from the world of integrable models, in [38j . 

Finally, we mention one further generalisation, first studied in detail in [32p39|: 
the 3-parameter family of T^T-symmetric Hamiltonians 

Hm,i,c. =P^- {ixf"' - aiixf^'-^ + 1(1 + l)/x2 (6.12) 

with boundary conditions on the Schrodinger-picture wavefunctions imposed in 
the same next-nearest-neighbour pair of Stokes sectors as before. The reality 
proof of [32j shows that the spectrum of TlM,a,i is 

• real for M > 1 and a < M + 1 + |2/ + 1| ; 

• positive for M > 1 and a < M + 1 - |2/ + 1| . 

This translates into the spectrum being entirely real in the regions below the 
dark dashed lines (red in colour) a = M + 1 it 2A on the plots in figure \T0\ 
where 2A = 2/ + 1. In fact, the full region of spectral reality is considerably more 
complicated than this initial result would suggest: the curved cusped lines on 
the plots (blue in colour) show the lines in the (2A, a) plane across which the 
number of complex eigenvalues changes, for various values of M. Further details 
will appear in ^40] . 

6.2 Physical consistency 

Given these surprising reality properties, a natural question arises: are VT- 
symmetric Schrodinger problems with entirely real spectra such as ()6.ip fully 
consistent quantum mechanical systems? Many recent papers [4H H2| B3l [HI H5] 
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Hi: I = -0.025 



Id: / = -0.0025 




Et: / = -0.0015 Hi: I ^ -0.001 



Figure 9: Reai eigenvaiues of — (ix)'^^^ + l{l+l)/x'^ as functions of M. 

(see also the review articles [MIIIT]) have addressed this question. Given a VT- 
symmetric Hamiltonian H in its unbroken VT phase, and denoting by {V'n} its 
set of eigenfunctions, up to a multiplicative rescaling we have 

VTMx) = M^) ■ (6.13) 
We can then define a VT inner product as 

{f,9)rr= [ dx[VTf{x)\g{x) , (6.14) 

where 7 is the contour in the complex plane along which the wavefunction is 
defined. However, the inner product (., .)-pT is not positive definite and it is easy 
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Figure 10: Lines in the (2A, a) plane across which complex pairs of eigenvalues appear, 
for various values of M . In each case, the spectrum is entirely real in the lower portion 
of the plot, and acquires complex eigenvalues as the curved lines are crossed. 

to check numerically that the orthogonality condition for the wavefunctions is 

(V'n,V'mW= (-l)"(5n,m , (n, m = 0, 1 , 2, . . . ) . (6.15) 
Similarly, indirect numerical support [HI [iHl [50l [51] leads to 

oo 

i'n{x)My) = {-irS{x - y) . (6.16) 

n=0 

Given the results (j6.15p and (j6.16p it is natural to introduce a new operator C 
such that 

V-^^ = CVTi^n = (-l)>n , (6.17) 

and define a CVT inner product 

{f,9)cvr= [ dx[CVTf{x)]g[x) . (6.18) 
77 
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One can check that a "PT-symmetric Hamiltonian H commutes with the CVT 
operator. Moreover, H is Hermitian with respect to the CVT inner product 
and it defines a unitary time evolution. The observables correspond to operators 
satisfying the following condition 

= CVTOCVT , (6.19) 

where stands for the transpose of O. A major difference between conventional 
quantum mechanics and PT-symmetric quantum mechanics is that the inner 
product (I6.18|) is generated dynamically by the Hamiltonian: in position space 
the linear operator C is given by 

oo 

c{x,y) = Y,Mx)My) ■ (6.20) 

n=0 

The parity operator V in position space is 

oo 

r{x, y) = J](-l)>„(x)V'„(-y) = 5{x + y) , (6.21) 

n=0 

and the time reversal operator T is just complex conjugation 

oo 

T{x, y) = j;(-l)"C(x)V^„(y) . (6.22) 

n=0 

Finally we note that further generalisations of VT symmetry have been considered 
- see for example the CVT quantum mechanical models discussed in [521 153j . 



7 The ODE/IM correspondence for the minimal mod- 
els 

The ODE/IM correspondence [23] connects the spectral properties of ordinary 
differential equations, including some of the "PT-symmetric problems discussed 
above, with objects which arise in the study of integrable models. The simplest 
example [23^ IM] links the following \d Schrodinger equation 

-^ + (x^--i?) + ^)^(x) = 0, (7.1) 

where M > and I is real, to the conformal field theories with c < 1 in the 
framework proposed by Bazhanov, Lukyanov and Zamolodchikov [551 156j . 

More precisely, equation (jT.lh encodes properties of the primary operator with 
conformal dimension 

(2; + l)i-4M^ 

^ ' ' 16(M + 1) ^ ' 



16 



belonging to a conformal field theory with central charge 



6M^ 

c(M) = 1 - i—— . (7.3) 

Give a pair of coprime integers p < q, the ground-state of the minimal model 
A^p^q is selected by setting 

M + l = ^, l + l = - (7.4) 
P 2 p 

in (j7.1|) . This corresponds to the central charge c(M) = Cpq = 1 — -^(q— p)^ , and 
lowest-possible conformal dimension 

MM,0 = -(l-(q-p)'). (7.5) 

pq 

For the Ising model A^3,4, equation ()7.5p gives h{^, — ^) = and the correspon- 
dence is with the identity operator, 1. For the Yang-Lee model M2.5, we have 
/i(|,— ^) = — ^ and the correspondence is with the relevant operator Now we 
observe that the singular term 1(1+1) /x"^ in (j7.ip with ^ = ^ — ^ can be eliminated 
by a change of variables 

X = , V(X, ^) = zP/4-l/2y(z, S) . (7.6) 

After rescaling z {2/p)'^^^z equation (17. ip finally becomes 

("^ + z^-\z^-P - E)^ y{z, E) = 0, (7.7) 

where 

E=n'-''^'E. (7.8) 



.2. 

The change of variable (|7.6p has replaced a singular potential defined on a multi- 
sheeted Riemann surface by a simple polynomial. Therefore any solution to the 
transformed equation (j7.7p is automatically single-valued around z = 0. 

Zero monodromy conditions have already played a central role in the ODE /IM 
correspondence [57^ i58j , and therefore it is natural to explore the constraints that 
this property imposes on the constant / in (j7.ip and (I7.2p [59] . 



7.1 Monodromy properties and the Kac table 

To see which other primary states have similarly-trivial monodromy, start from 
()7.ip with generic / > —1/2 and perform the transformation (17. 6p . The result is 

f-^ + + z'^-'iz^-^ - E)) y{z, E, [) = , (7.9) 
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where 2{l + i) = p(/ + i). The presence of the Fuchsian singularity at z = 0, 
compared to (j7.7p , means that the zero monodromy condition of a generic solution 
is no longer guaranteed. In general, equation (17.90 admits a pair of solutions 



oo 



n=0 n=0 

where Ai = / + 1 > A2 = — ^ are the two roots of the indicial equation A(A — 1) — 
[([+ 1) = and 

Xj{e^-'z) = e^^'^^Xj{z) , J = 1,2. (7.11) 
Therefore, the general solution is a linear combination 

yiz,Ej) = axiiz)+TX2iz) (7.12) 

and we shall demand that the monodromy of y{z) around z = is projectively 
trivial: y{e^'^^z) oc y{z). We shall see that this condition imposes 

i) 2/ + 1 is a positive integer; 

ii) The allowed values of 2/ + 1 form the set of holes of the infinite sequence 

pr + qs , r,s = 0,1,2,3... . (7.13) 
We call the set of integers (17.13j) 'representable', and denote them by Mpq. 
As a consequence of i) and ii) , it turns out that the allowed conformal dimensions 

HM,l) = ^l±^l^^^ (7.14) 
4pq 

exactly reproduce the set of conformal weights of the primary operators in the 
Kac table of M p,q . Figures [TT] and [12] illustrate the situation for the Ising and 
the Yang-Lee models. 

^ ^ 
16 2 

• oo**o«««« 

0123456789 



Figure 11: The holes in the infinite sequence of integers for the critical Ising model 
A^3^4. The holes are at 1, 2 and 5. 

To establish these claims, we first notice that the requirement that the general 
solution y{z) be projectively trivial means that Xiiz) and X2{z) must have the 
same monodromy. This implies 

Ai - A2 = 2[+ 1 G N, (7.15) 
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5 

• o«o**«««« 

0123456789 

Figure 12: Holes for the Yang-Lee model 7^2,5 , at 1 and 3. 

and therefore restricts / to be an integer or half integer. However, in such a cir- 
cumstance, while Xi{^) keeps its power series expansion, X2iz) generally acquires 
a logarithmic contribution 

^ oo 

X2{z) = Dxi{z) log(z) + - ^ dn^" , (7.16) 

^ n=0 

and unless D = 0, this will spoil the projectively trivial monodromy of y{z). The 
log term is absent if and only if the recursion relation for the dnS with D = 

n{n- 21 -l)dn = dn-q - E dn^p (7.17) 

with the initial conditions do = 1, dm<o = admits a solution. 

Consider first the case 21 + 1 ^ Mpq. Starting from the given initial conditions, 
the recursion relation generates a solution of the form 

^ oo 

X2{z) = -^dnZ^ (7.18) 

^ n=0 

where the only nonzero d^s are those for which the label n lies in the set Mpq. 
Given that 21 + 1 ^ IRpq, for these values of n the factor n{n — 21 — 1) on the 
LHS of the recursion relation is never zero, and the procedure is well-defined. If 
instead 21 + 1 £ Kpq, then the recursion equation taken at n = 2/ -|- 1 yields the 
additional condition 

^Wp- Wq=0' (7-19) 

which is inconsistent for generic E, and so the log term is required. If the coprime 
integers p and q are larger than 1 then the number of 'nonrepresentable' integers 
Npq with Z"*" = Mpq U Npq is in fact |Npq| = ^(p— l)(q— 1), a result which goes 
back to Sylvester [60j and matches the number of primary fields in the Kac table. 

8 Applications to condensed-matter physics 

There are many interesting potential physical applications of PT-symmetry, rang- 
ing from optics to the more speculative proposals related to conformal gravity |61| . 
and Higgs phenomena in the Standard Model of particle physics [62j. A partial 
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list of examples can be found in the review [37| on Pseudo-Hermitian quantum 
mechanics. An 'experimental' means of distinguishing between a (Dirac) Her- 
mitian and a "PT-symmetric Hamiltonian through the quantum brachistochrone 
is an interesting recent development [63^ [M| . In section [3] a general relationship 
between the Yang-Lee model and non-intersecting branched-polymers was dis- 
cussed. Here we would like to report on a simple application of "PT-symmetric 
quantum mechanics to the ferromagnetic Ising system and briefly mention a sur- 
prising correspondence between the equation (jT.ip and experiments in cold atom 
physics. 

Recently Lamacraft and Fendley [65] were able to connect the study of uni- 
versal amplitude ratios (see also [66]) in the quantum spin-chain version of the 
Ising model to a simple PT-symmetric Schrodinger problem. The Hamiltonian 
of the quantum spin-chain Ising model is 

L 

n\a, h] = - Yj^at + af af+ J , (ai+i = af ) , (8.1) 



where and o"^ are Pauli matrices. The model (j8.ip is critical at h = he = 1 
and in its ordered and disordered phases for h < 1 and h > 1 respectively. 
The magnetization operator M = X^i'^l^/^ does not commute with 'H. Setting 
h = he and considering the continuum limit versions of the operators Ti and 
M = Jq dx a{x), we can define the magnetisation distribution in the ground- 
state |0) as 

P{m) = mm - M)|0) = / ^e-'^^xW , (8.2) 

where 

X(A) = (0|e^^^|0) (8.3) 
is the generating function of the moments of the distribution P{m) 

d 



(0|M |0) = (m2") = / m2"P(m) dm = (-1 

J —oo 



rx(A) 



(8.4) 



It is a direct consequence of the Yang and Lee results |17[ [T5] that the function 
x(A) is an entire function of A, and a consequence of the Hadamard factorisation 
theorem that 



n=\) ^ 



(8.5) 



Finally x can be related to the spectral determinant T |37j of the following PT- 
symmetric Hamiltonian (cf. (j6.1ip ) 

H^^_,=p'-{^xr-^^. (8.6) 
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The relation is simply 

X(A) = ^T{E) , X^^E . (8.7) 

The second application, proposed by Gritsev et al. in ji67j , is related to inter- 
ference experiments in cold atom physics. In a typical interference experiment 
[67l EHl Eh] , a pair of parallel Id condensates of length L is created using standard 
cooling techniques and an highly anisotropic double well radio-frequency-induced 
micro-trap. After the condensates are released from the trap, interference fringes 
are observed using a laser beam. If the condensate has the form of two coaxial 
rings [69], the associated 'full distribution function' in the zero temperature limit 
can be mapped, through the ODE/IM correspondence, to the spectral determi- 
nant of the Schrodinger equation (j7.ip . 

9 Conclusions 

In this short review we have discussed non-Hermitian quantum mechanics, in- 
tegrable quantum field theory, links between them via an ODE/IM correspon- 
dence and some applications to condensed-matter physics. The material is by no 
means exhaustive but rather refiects the authors' experience and interests. Apart 
from a small number of exceptions, most of the papers on T^T-symmetry concern 
one-dimensional Schrodinger problems in the complex domain; many fundamen- 
tal issues related to their full consistency as quantum mechanical models have 
been successfully addressed. Since the early 1990's two-dimensional relativistic 
non-unitary quantum field theory, particularly the 2d scaling Yang-Lee model, 
has been extensively studied by the integrable model community. Considering 
the relatively-restricted physical relevance usually given to these models, a large 
quantity of exact results have been produced over the years. However, none of the 
fundamental interpretation issues, of great concern in the T^T-symmetry commu- 
nity, have been completely clarified. We think this in an interesting open problem 
which deserves attention. 
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